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ABSTRACT 



> 

o 
o 

\^ ■ When a self-gravitating spherical gas cloud collapses or accretes onto a central 

^^ ■ mass, the inner region of the cloud develops a density profile p oc r^'^' ^ and the velocity 

Q^^ . approaches free-fall. We show that in this region, nonspherical perturbations grow with 

(— I , decreasing radius. In the linear regime, the tangential velocity perturbation increases as 

S^l r~^, while the Lagrangian density perturbation, Ap/p, grows as r~^'^. Faster growth 

O ' occurs if the central collapsed object maintains a finite multiple moment, in which case 

Ap/p increases as r~ , where / specifies the angular degree of the perturbation. These 
C^ ■ scaling relations are different from those obtained for the collapse of a homogeneous 

cloud. Our numerical calculations indicate that nonspherical perturbations are damped 
^ . in the subsonic region, and that they grow and approach the asymptotic scalings in the 

supersonic region. The implications of our results to asymmetric supernova collapse 
and to black hole accretion are briefly discussed. 



Subject headings: hydrodynamics - accretion - supernovae 



1. Introduction 

Gravitational instability is responsible for a wide range of structures in the observable 
universe. Much effort has been devoted to structure formation in cosmological models, and it is 
well-established that the Hubble expansion slows the growth of perturbations (e.g., Peebles 1980). 
In this paper we are interested in the behavior of nonspherical perturbations in the background of 
spherical collapse or accretion. 



Hunter (1962) studied the collapse of a homogeneous, pressureless, dust cloud, and showed 
that in the linear regime, perturbations of arbitrary shape and scale grow asymptotically as 
5p/ p oc (to — t)~^ oc i?~'^'^, where t = t^ denotes the time of complete collapse, and R is the 
cloud radius. It was thought that this instability might be responsible for the fragmentation of 
collapsing protostellar clouds. Lin, Mestel and Shu (1965) studied the collapse of a homogeneous 
ellipsoidal cloud of dust, and showed that its ellipticity increases until a sheet or pancake forms. 
Homogeneity plays a key role in both examples. However, the presence of a slight initial central 
concentration significantly alters the evolution of a cloud. Since the dynamical time, (Gp)~^'^, 
is shortest in the central region, a cuspy density and velocity profile develops. As a result, 
one may question the applicability of the Hunter and Lin-Mestel-Shu instabilities to realistic 
astrophysical situations. Goodman &; Binney (1983) have already commented on problems with 
the Lin-Mestel-Shu instability for inhomogeneous clouds. In this paper, we show that perturbation 
growth is slowed by the central mass concentration in an inhomogeneous collapse. 

A related problem concerns the stability of spherical accretion flow. Bondi (1952) found a 
class of solutions describing steady-state accretion onto a compact object from a homogeneous 
medium. There is a unique transonic "critical" flow, for which the mass flux is maximum; the 
other subcritical solutions describe subsonic flows. Bondi speculated that nature would prefer the 
critical flow (see Shu 1992 for a discussion) and underlined the importance of a linear stability 
analysis. After several attempts by a number of authors, the correct analysis was achieved by 
Garlick (1979) and Moncrief (1980), who showed that both the critical and subcritical flows are 
globally stable. However, this does not mean that perturbations can not grow spatially. Indeed, 
we show in this paper that nonspherical perturbations carried by fluid elements are amplifled in 
the supersonic region of the critical flow (see also Goldreich, Lai & Sahrling 1996; Kovalenko & 
Eremin 1998). 

The present study originated from our attempts to understand the origin of asymmetric 
supernova (Goldreich, Lai & Sahrling 1996; hereafter GLS). A large body of evidence suggests 
that Type II supernovae are globally asymmetric, and that neutron stars receive kick velocities 
of order a few hundred to a thousand km s~^ at birth (see, e.g., GLS, Cordes & Chernoff 1998, 
and references therein). The origin of the kick is unknown. A class of mechanisms relies on local 
hydrodynamical instabilities in the collapsed stellar core (e.g.. Burrows et al. 1995; Janka & 
Miiller 1994, 1996; Herant et al. 1994) that lead to asymmetric matter ejection and/or asymmetric 
neutrino emission; but numerical simulations indicate these instabilities are not adequate to 
account for kick velocities ^ 100 km s"^ (Burrows & Hayes 1996; Janka 1998). Global asymmetric 
perturbations of presupernova cores may be required to produce the observed kicks (GLS; Burrows 
& Hayes 1996). GLS suggested that overstable g-modes driven by shell nuclear burning might 
provide seed perturbations which could be amplifled during core collapse.^ 



^See Lai & Qian (1998) and Arras & Lai (1999) for discussion/review on alternative mechanisms which rely on 
asymmetric neutrino transport induced by strong magnetic fields. 



Hints regarding perturbation growth during collapse are obtained by considering the collapse 
of finite-mass fluid shells onto central point masses. In the context of core-collapse supernova, 
one might imagine that the spherical shell mimics the outer supersonic region, while the central 
mass represents the homologous core. A large number of growing modes can be identified (GLS; 
Appendix B). In particular, the dominant dipole mode grows according to AR{6, t)/R{t) oc {tQ — ty, 
with s < 0, where R is the mean radius of the shell and Ai? is the perturbation. For a core-shell 
mass of order unity, the power-law index s ~ — 1, corresponding to AR/R oc R^^'"^ . 

The thin-shell model is too simplistic to be applicable to real presupernova collapse. One 
approach is to determine the stability of the self-similar collapse solution (Goldreich &; Weber 
1980; Yahil 1983). An analysis by Goldreich & Weber (1980) shows that the inner homologous core 
is stable against nonradial perturbations. This is not surprising given the significant role played by 
pressure in the subsonic collapse. Pressure is less important in the supersonically collapsing region, 
making it more susceptible to large scale instability. A stability analysis of Yahil's self-similar 
solution, which extends the Goldreich- Weber solution to include a supersonically collapsing outer 
core, does not reveal any unstable global mode before the proto-neutron star forms (Lai 2000). 
However, it is perhaps more illuminating to treat the initial-value problem, and determine how 
initial perturbations evolve as the collapse proceeds. We carry out such an analysis in this paper. 

The rest of the paper is organized as follows. The basic perturbation equations are 
summarized in §2, and in §3 we derive the asymptotic scaling relations for the perturbations in 
the regime where the collapse/accretion flow is supersonic. In §4 we present a numerical study of 
the evolution of perturbations during collapse; this numerical study not only confirms the analytic 
asymptotic relations, but also explores the regime where pressure is important. We discuss the 
implications of our results in §5. 



2. Basic Equations 

We consider barotropic fluid obeying the equation of state p = Kp'' , where 7 is the adiabatic 
index and K \s a, constant. The unperturbed flow is spherically symmetric, with velocity in the 
radial direction. The Eulerian perturbations of density, /?, and velocity, v, can be decomposed into 
different angular modes, each of which has the form 

5p{v,t) = 6p{r,t)YiUe,(t>), (1) 

,5v(r,t) = 6vr{r,t)YUOA)r + 5v^ir,t)V^YU0A)+V±x[5v,ot{r,t)YUe,(l))f], (2) 



where 



Vx^4 + ^|-, (3) 

ov sm 00 



and f, 9, (p are unit vectors in spherical coordinates. The perturbations of pressure, p, and 
gravitational potential, ip, have the same angular dependence as 6p{r,t). The perturbed mass 



continuity equation reads 

fi + (V.v)Jp+^(,.V*)'-'('+l)^=«. (4) 

where p, v = vr denote the unperturbed spherical flow variables, prime stands for d/dr, and 
d/dt = d/dt + V • V is the total time derivative. The perturbed radial Euler equation can be 
written as 

''''^+v'5v. = -M-m'. (5) 



dt \ p 

and the tangential Euler equation reduces to 



''^^ + -8v^ = -- {'A - -6i., (6) 



dt r r \ p J r 

|(r5^rot) = 0. (7) 



The perturbed Poisson equation is 

J_ [r\6i^y] ' - ^^^X^h i^ = A7rG6p. (8) 

The vorticity perturbation reads 

V x6^r = f X V_Ll/m H SvrotYim r + -{rdvrot)' ^±Yim, (9) 

/va fyt tyt 

where 

5vt = Svr — Su' , (10) 

with 

5u{r,t) = r5v±{r,t). (11) 

Thus (5frot and 6vt are related to the vorticity of the perturbed flow. Equation (|6|) is transformed 

to 

dSu 6p 

-— = 5ip, (12) 

at p 



with the aid of equation (11). Combining equations (g) and ([l^) , we flnd 

-—5vt = —v'5vt (13) 

dt 

Equations (0) and (13) express the conservation of circulation in a barotropic fluid.0 The 
former also reflects the conservation of angular momentum; following a fluid element, 6vmt oc 1/r. 
Since, as we shall prove shortly, v oc r~^'^ as r ^ 0, equation ( [T^ ) implies the asymptotic relation 



^For homogeneous collapse or expansion, as in an expanding universe, we have v = {R/R)r, where R is the scale 
factor. Equation (|l3) then becomes the familiar d{R5vT)/dt — 0. 



5vt oc r^'^. We shall focus on irrotational flows from here on. We neglect 6vj-ot because it is 
decoupled from the density perturbation, and 6vt because it decays inward. The continuity 
equation (^) for irrotational flow simplifies to 

d5p 
~dt 



+ (V • w)5p + ^ {r^pSu')' - 1(1 + 1)^ = 0. 



Note that the Poisson equation (S) has the following integral solution: 



<5V(r,i) 



AttG 
21 + 1 



1 



r-^+l 



Qiir,t) + r'Siir,t) 



where 



Qi{r,t)= / x^+^5p{x,t)dx, Si{r,t) 
Jo 



X Sp{x, t) dx. 



(14) 



(15) 



(16) 



Equations (p^, (14), and (p^ ) determine the perturbed flow. This particular form of the 
perturbation equations is convenient for implementation as a Lagrangian numerical code (see §4). 



3. Asymptotic Analysis 

Consider a cloud in hydrostatic equilibrium with an initial density profile that decreases 
outward. As its pressure is depleted, the cloud starts to collapse. Since the dynamical time, 
(Gp)~^'^, decreases outward, cuspy density and velocity profiles will be established after the core 
has collapsed. By contrast, a uniform density dust cloud collapses homologously, and remains 
uniform as the collapse proceeds. We now study the behavior of the flow and its perturbations in 
the asymptotic regime where the gas pressure is negligible compared to gravity. 



3.1. Unperturbed Spherical Flow 

Consider how the velocity, Vm^ and density, pm, of a fluid shell with enclosed mass m change 
as the shell collapses from its initial radius Vmo to a smaller radius r^- The pressure is negligible 
in the supersonic region, so for r^ <C rmo, we have 



f2Gm\ 

V Tm / 

Then from the continuity equation, we obtain 



oc — r, 



^1/2 



Pm OC r, 



-3/2 



(17) 



(18) 



Note that these relations do not require the mass of the collapsed object to be fixed. Indeed, with 
mc{t) the mass of the collapsed core at time t, we have 

"2Gmc(t)T^/^ 



v{r,t) 



oc r ' , p{f,t) ~ 



mc{t) 



oc r 



-3/2 



(19) 



where rhc is the mass accretion rate onto the core.^ When the accretion time, rric/rhc, is 
much longer than the dynamical time of the flow, equation ( [l9| ) describes the inner region of a 
steady-state Bondi flow. In a dynamical collapse, when the accreted mass becomes much larger 
than the original core mass, dimensional analysis implies rndt) = i^^'^G~'^'^~^"^t^~^'^mc, where 
fhc is a dimensionless number, and t is measured from the moment when the center collapses. 
For 7 = 1, this reduces to the familiar rric = rhc{cl/G)t (see Shu 1977). Equation (|T9|) describes 
the central region of Shu's expansion- wave solution (Shu 1977), and the post-collapse extension 
of the Larson-Penston solution (Larson 1969; Penston 1969; Hunter 1977) in the context of star 
formation, as well as Yahil's post-collapse solution in the context of core-collapse supernova (Yahil 
1983). 

Note that the above asymptotic scaling solution assumes supersonic flow for r ^ 0, i.e., 
V ^ Cs oc p(f~^)''^ oc r"^^'*'"^-''^. This requires 7 < 5/3. The special case of 7 = 5/3 is considered 
in Appendix A. 



3.2. Perturbations 

We investigate asymptotic power-law solutions to equations (|T^), (p^) and (15). Let 6p ex r° 
and 6u oc r . The Poisson equation has a general solution of the form 6ip ~ Qc/'^ + r'^Sp, where 
the first term arises from a central multipole moment, Qc, and the second term is due to the 
density perturbation outside the the central core. In most astrophysical situations, Qc is zero or 
close to zero; possible exceptions are discussed separately in §3.3. For example, in accretion onto 
a star, the supersonic flow may be stopped by a standing shock near the stellar surface. Inside the 
shock, any inhomogeneity carried in by the gas will be smeared out on a local dynamical timescale. 
In accretion onto a Schwarzschild black hole, the event horizon defines the inner boundary of the 
supersonic flow, and the no-hair theorem ensures that mass multipole moments are not retained 
by the black hole. Thus we have 6ip oc r^^"^. 

In the asymptotic regime, d/dt — > v{d/dr), so equations (12) and (|l^) reduce to: 

^ + -fKp-y-^6p + # = 0, (20) 

^ = 0. (21) 



6-l|-i(i + l) 



3\ v6p 
"+2)V^ + 

Equation ( |2lD implies b = a + 2 for a 7^ —3/2. Equation ( ^0[) has the scaling form 
0{br^-^/^) + C'(i^r"-3{7-2)/2) _^ 0{r''+^) = 0, from which we see immediately that b = and 
a = — 2 (for 7 < 5/3). To obtain the scaling behavior for Sv^ = Su' , we need a higher order 



■^If the central region of a pressureless cloud is non-singular to begin with, then at the moment when the center 
reaches infinite density, the central density and velocity profiles are p oc r'^^''^ and v oc r'^'"' (Penston 1969). 
However, after the core has formed, the profiles are given by equation (|l9|). 



correction for 5u. Let 5u = 5uq + 5uir^^ , where 5uo and 5ui are constants independent of r. For 
K = (the pressureless case) or 7 < 2/3 so that 5p/ p < dip asymptotically, equation (|l^ gives 
0{biv6uir^^-'^) ~ 0{r''+^). We then have bi = a + 7/2 = 3/2. For iT / and 7 > 2/3, equation 
(|l§) reduces to 0{bivduir''^~'^) + C'(r"-3(')'-2)/2) = q, which gives 61 = (5 - 37)/2. To summarize, 
the asymptotic scaling relations for the perturbations are: 



rv 
where 



5u ~ 6uQ + 6uir''\ (22) 

5vr oc r^i-\ (5t;^ocr-\ (23) 

5p ~ -2/(/ + l)^ocr-2, (24) 



^ r3/2 forK = 0or7<2/3, 

^ 1 (5 - 37)/2 for K / and 7 > 2/3. ^ ' 

Note that the above results apply for 7 < 5/3. The special case of 7 = 5/3 Bondi accretion is 
discussed in Appendix A. 



3.3. Physical Interpretation 



Equations (22)-(^) describe the Eulerian perturbations. To derive scaling relations for the 
Lagrangian displacement, ^(r,t) = S,rir,t)Yimf + S,±{i"^t)^±Yimi and for the Lagrangian density 
perturbation, Ap, we note that Av = (5v + (^ • V)v = d^/dt = (d^/dt) + (v • V)^, which yields 

SVr = ^+<-^'Cr, (26) 

Sv± = ^ + <l--e±. (27) 

ot r 



In the asymptotic regime, we find using equation ( p3D that 

The Lagrangian density and velocity perturbations are 

^ ~ ^ oc r-y^ ^ oc r''-V2. (29) 

p V V 



The scaling relations (|28| ) and (^9D can also be derived directly from Lagrangian perturbation 
theory. 



Equations (^2|)-(p4|) and ( |29[) are the main results of this paper.Q| They can be understood 
from the following simple consideration: The time that a fluid element spends near radius r, of 



This result was discussed without derivation in GLS. Similar results were also obtained by Kovalenko & Eremin 
(1998) in the context of Bondi accretion. 



order tr = dt/dlnr oc r^'^, decreases rapidly as r decreases. Unless the specific torque, [Sil^ + Sp/p], 
increases inward sufficiently rapidly to compensate for the decreased time, the specific angular 
momentum, Su, of the fluid element will be independent of r at small radii. Thus Av± oc r~^. 
This leads to a tangential displacement ^_l ~ trSv± oc r^'^, which induces a density perturbation 
Ap/p ~ ^±/r oc r"^/^. 



Equation (p9[) indicates that the fractional density and tangential velocity perturbations grow 
in supersonic collapse/accretion. The radial velocity perturbation, however, can be affected by 
pressure even in the regime where the pressure has negligible effect on the unperturbed flow. We 
see that Avr grows with decreasing r when 7 > 1, and Av^/v grows only when 7 > 4/3 (but recall 
that the scaling relations apply only for 7 < 5/3). 



3.4. Special Cases: Possibility of Faster Growth 

Now consider a hypothetical situation in which the inner boundary of the flow is a "sticky 
sphere": Once a fluid element enters the sphere, it gets stuck on the spot where it enters. In this 
case, a finite multipole moment, Qc, will accumulate at the center (see eq. fl^). In the region 
where Qc dominates the gravitational perturbation, we have 5ip ~ Qc/r^^^ ^ r'^5p. It is easy 
to show from equations (pOD-(^Il) that the perturbations have the following scaling behavior (for 
7 < 5/3): 

5u ~ —^-54) o^ r^''^~\ (30) 

21 — Iv 

5p ^ -^^ocr-'-=^/^ (31) 

47rG \ Qc __ „_;-! 



'* = - (27TIJ ?Tr « -■"■■■ (32) 

These scalings depend on /. Even for / = 1, the growth is faster than the case discussed in §3.2 
and §3.3. These scalings can be understood as follows: The central multipole moment exerts a 
torque oc r^'^^ and a radial force oc r^'^^ on a fluid element. The angular momentum grows 
as trSip oc r +^1'^ ^ and the radial velocity perturbation grows as 6vr oc r~ /^. Thus we have 
5vr/v ~ 6v±/v oc r~K The Lagrangian displacement is Cr ~ '^± ~ trSv oc r~'+^, and the resulting 
density perturbation is Ap/p ~ ^/r oc r~K 

As we have argued in §3.2, it is unlikely the scaling relations derived in this section will apply 
in general situations since it is hard to imagine that the central core can retain its multipole 
moments for a time much longer than the dynamical time of the flow at the inner boundary. There 
are two possible exceptions: 

(i) For I = 1 modes: The core has an extra degree of freedom, i.e., it can have linear motion in 
response to flow perturbations. This gives rise to a net central dipole moment and the possibility 
of faster growth of dipolar perturbations. Let rric be the core mass and Zc be the position (along 



the z-axis) of its center-of-mass relative to the the origin of coordinates. The displaced core 
induces a potential perturbation 6ip ~ —{GrricZc/r'^) cos 6. From equations (|30|)-(p^, the flow 
perturbations are given by 

5p 3Zc Svr Zc 5v± Zc . , , 

— ~ cos 61, ~ — cos 6^, ~ — smb', (33) 

p 2r V 2r v r 

where we have used v ~ — (2Gmc/r)^'^. However, the asymptotic perturbations given by equation 
( p^ ) simply represent a spherical flow centered at Zc- Thus it is not surprising that the sum of the 
gravitational force exerted on the core plus the rate at which momentum flows across a surface 
surrounding it sum to zero. 

(ii) For / = 2 modes: Suppose the central core consists of a rotating star, or a star with 
a massive circumstellar disk. This could result from an early phase of accretion/collapse with 
significant angular momentum. Subsequent spherical accretion, with no net angular momentum, 
will be affected by the central quadrupole. 



4. Numerical Calculations of Linear Perturbation Growth 

The asymptotic scaling relations derived in §3 apply only in the supersonic regime. To 
determine the behavior of perturbations under general conditions, we numerically follow the 
collapse of a self-gravitating cloud and evolve the nonspherical perturbation carried by each 
fluid element. Because of the large disparity in the timescales involved in the central region 
and the outer region, it is essential for the code to have a wide dynamical range. Previous 
multidimensional simulations of Bondi accretion (e.g., Ruffert 1994) did not achieve high enough 
resolution in the central region to reveal the growth of perturbation. We restrict our calculations 
to the linear regime. Thus perturbations associated with different Y/^ evolve independently, and 
our calculations involve only one spatial dimension. 

We have constructed a one-dimensional Lagrangian finite-difference code. The unperturbed 
flow variables {r,v,p) are followed with a standard scheme (Bowers &: Wilson 1991), and the 
Eulerian perturbations {6p,5u,6ip) are evolved using equations (|l^), (|14|) and (p^. The flow is 
covered by a uniform mass grid. The quantities r, v, 6u, dif) are zone-edge-centered, while p, 5p 
are zone-centered. A staggered leapfrog integration scheme is adopted to ensure second-order 
accuracy in time. 



4.1. Pressureless Collapse 

To calibrate our code and check the asymptotic scaling relations of §3, we study the collapse 
of a centrally concentrated dust cloud. Figure 1 shows an example of such a collapse calculation. 
The cloud, of total mass m = 1 and radius r = 1, is initially at rest with density proflle p oc r~^; 



- 10- 

i.e. the radius of a mass shell with enclosed mass m is r^ = 'm}''^- We initialize an / = 2 
perturbation with 6p/ p = 1 in arbitrary units and 5u = 0. The Poisson equation is solved to 
give the potential perturbation, 5^. We impose an inner boundary at r^ = 10~^: Once a mass 
shell enters this boundary, it is removed from the simulation domain, and its perturbations are 
immediately smeared out; i.e., the central quadrupole moment, Qc, is maintained at zero. The 
profiles p oc r^^' ^ and v oc r^^' ^ for the unperturbed flow are established near the center as the 
collapse proceeds. Figure 1 depicts the evolution of the perturbation carried by three different 
mass shells. We see that as the mass shells collapse to small radii, the analytic asymptotic scalings 
5u — > constant and 5p/ p oc r~^' ^ are achieved. Calculations with other initial conditions confirm 
that these scalings are generic features of perturbation growth in the absence of central multipole 
moments. 

We have also studied the case where the central multipole moment, Qo is nonzero, and 
confirmed the steeper scalings derived in §3.4. 



4.2. Collapse with Finite Pressure 

Next we study the collapse of clouds having finite pressure. For definiteness, we choose the 
initial cloud to be a 7 = 4/3 spherical polytrope in hydrostatic equilibrium. The collapse is 
initiated by reducing 7 to 1.3, and by reducing K by 10%. Either one of these reductions alone is 
adequate to induce the collapse. This mimics Type II supernova collapse, where a white dwarf core 
of a massive star collapses to a neutron star. However, to focus on the growth of perturbations 
during the collapse, we do not include a shock in our calculation. We define an inner boundary 
at Tc = 0.005; the original cloud radius is r = 1 and its mass m = \. When the central density 
becomes greater than 10^, corresponding to a few times nuclear density if the initial cloud has 
mass and radius typical of a Chandrasekhar mass white dwarf, we cut out the flow inside r^ from 
the computational domain. This enables us to follow the collapse and accretion of the rest of the 
cloud. 

Figure 2 shows the unperturbed density and velocity of several different mass shells as 
functions of their Lagrangian radii. The inner region collapses homologously; the Mach numbers 
of individual shells remain below unity outside Vc- The outer region of the flow goes through a 
transonic point, and eventually attains the free-fall asymptotics, with v oc r^^'^ and p oc r^^'^. 
Inspecting the flow profiles at different time slices (not shown), we confirm that our numerical 
results agree with the self-similar solution (Yahil 1983) in the regime where it applies. 

Figure 3 and Figure 4 show two examples of the evolution oi I = 2 perturbations during the 
collapse depicted in Fig. 2. In Fig. 3, the initial perturbation is chosen to be 6p/ p = 1 and 5u = 0, 
while in Fig. 4, the initial perturbation corresponds to the eigenfunctions of the first g-mode of a 
7 = 4/3 polytrope, with adiabatic index 71 = 5/3. This value of 71 is used only for setting up the 
initial perturbations; after the collapse starts, the adiabatic index is set to 7 = 1.3. We see that 
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the perturbations carried by the inner mass shells {m = 0.2, 0.6), which never become supersonic, 
vary in an oscillatory manner with no increase in amplitude. This is consistent with the result of 
Goldreich & Weber (1980) that the homologous inner core of a collapsing 7 = 4/3 polytrope is 
stable against non-radial perturbations. However, the outer region of the cloud {m ^ 0.85) attains 
a high Mach number and eventually approaches free-fall. We see from Figs. 3 and 4 that the 
density perturbation grows in this outer region, and that the asymptotic scaling relations derived 
in S3. 2 are recovered. 



5. Discussion 

Nonspherical perturbations are amplified during the supersonic collapse or accretion of a 
centrally concentrated gas cloud. We have derived asymptotic scaling relations for their growth. 
These general results have implications for several different astrophysical problems which we now 
discuss. 

In spherical accretion onto a black hole, we expect the radiative efficiency to increase as a 
result of nonradial perturbations in the flow. In the asymptotic regime, the tangential velocity 
perturbation scales as 6vj_ oc r~^, and the corresponding Mach number scales as 6v±/cs oc A^'^'~'^'''^. 
For 7 < 5/3, the Mach number grows faster than r"^'^. Similarly, the Mach number associated 
with the radial velocity perturbation (see eq. [23| for 7 > 2/3) scales as 6vr/cs oc r~^^'^~^'''^. 



One might expect the formation of shocks which lead to thermalization of the flow and high 
radiative efficiency (see Chang & Ostriker 1985 for previous discussion on the formation of shocks 
in spherical flows). This point has also been noted recently by Kovalenko & Eremin (1998), who 
derived similar scaling relations for Bondi accretion. 

In the context of core-collapse supernovae, our results indicate that perturbations in the 
homologous inner core do not grow, but those in the outer core, involving ~ 15% of the core mass, 
and in the envelope are amplified. Since 6p/p scales as r~^'^, we expect that the amplification 
factor is at most 10 for r decreasing from 1500 km to 15 km. It is possible that dipole perturbations 
obey the r~^ scaling (see §3.4), in which case the amplification factor could be larger. Interestingly, 
if overstable g-modes driven by shell nuclear burning are responsible for the seed of presupernova 
perturbations (see GLS), it is exactly at the outer core where the perturbation amplitude is 
expected to be the largest. The asymmetric density perturbation may lead to asymmetric shock 
propagation and breakout, which then give rise to asymmetry in the explosion and a kick to the 
neutron star (e.g.. Burrows & Hayes 1996). 

Finally, we note that our analysis has neglected rotation (i.e., the flow does not have a net 
angular momentum). In the perturbative regime, rotation (around the z-axis) is represented by 
the last term of equation (^) with m = 0, i.e., (Jvi-ot = —^Viot{f,t){dYiQ/d6)(j). As shown in §2, 
this rotational perturbation is decoupled from the density perturbation. Therefore, provided the 
rotational velocity is small in comparison to the radial velocity, i.e., \6vt\ <C \v\, we expect our 
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scaling relations for the growth of perturbations to be valid. 

A major portion of this research was done between 1995 and 1997, when D.L. was a postdoc 
in theoretical astrophysics at Caltech; support from a Richard C. Tolman fellowship is gratefully 
acknowledged. This work is supported in part by NASA grant NAG 5-8356, by a research 
fellowship from the Alfred P. Sloan foundation, and by NSF grant 94-14232. 



A. Perturbation of Bondi Accretion for 7 = 5/3 

The Bondi solution for 7 = 5/3 is special because the flow remains subsonic for r > 0. We 
can think of the sonic point as being located at r = (in Newtonian theory). For r <^ GM/c^, 
where M is the central mass and Cqo is the sound speed at infinity, the flow velocity and density 
are given by 



T^ , (Al) 



\ 2r J 
'GM\ 

Although these have the same scalings as equations ([l7|) and (|l^ ) , the asymptotic behavior of the 



~ " [iTir] r-'"- (A2) 



perturbations is quite different. Indeed, for 7 = 5/3, eqs. ( pOl) and (21) yield 



5pocr'-^/^ (57iocr'+3/^ ^V^ oc r'. (A3) 

Thus all perturbations decrease as r decreases in keeping with the subsonic nature of the 
unperturbed flow. If a central multipole moment is present (see §3.4), the following asymptotics 
become dominant: 

Jpocr^'-^/^ <5ucxr-'+^/2^ <5V oc r^'^^ (A4) 



B. Perturbations in Collapsing Spherical Shells 

Consider a spherical fluid shell of mass Mg falling from infinity onto a central point mass M^- 
The shell radius, Ro{t), evolves in time according to 

d'^Ro GM . ., ., Ms 

_o=__ w,th M = M. + -, (Bl) 



which gives 

Ro{t) = [^) i-tr', (B2) 



f9GM\^/\ 2/3 
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where we have set t = at the point of complete collapse. The surface density and radial velocity 
are given by 



So(t) 



^' -' ^^-4/3 T/.m = ^ ^ r-f^-1/3 



^-nRl 



« (-t)- 



y^it) 



dt 



cc i-ty 



(B3) 



B.l. Perturbation Equations 

The dynamical variables for a perturbed shell are its surface density S(^,t), radius R{6,t), 
radial velocity Vr{0, t) = R, where dot indicates d/dt, and tangential velocity V^(0, t). We assume 
spherical coordinates and axisymmetry so that there is no (/)-dependence. The core mass, M^, 
is free to move. We use Zc{t) to denote its displacement from the coordinate origin. Note that 
S, Vr, V_L are rigorously defined from the three-dimensional fluid variables via 



R^Jr- 



pr dr, 



Vr = 



1 



i^^ 



pr Vrdr, 



V, = 



1 /""^ 2 , 



where the integration runs through the thickness of the shell. Using these definitions and the 
standard hydrodynamical equations, we derive the continuity and Euler equations for the shell: 



(B4) 



as _ 2R 
'dt ~~"r" 



dVr 



R 



SVx-V^, 



RV± 



1 r 



+ 



dr / R^ 



(Vx'I>)^+ + (Vx<&)h_ 



dt R 2 

lat the nonspherical perturbation is small. Note that here 
e {d / 89) + {4>/ sin e){d/d^)] = {l/R)e{d/de). The gravitational 



(B5) 
(B6) 

(B7) 



where we have assumed t 

Vx = (l/i?)Vx = {l/R) 

potential <J> = <I>c + ^s includes contributions from the core, <l>c) ^^ given by 

GMr „.. ^^ Zi ( 47r \i/2 



^c(r,t) 



ZcZ\ 



-GM,^:.^' 



'+1 V2/ + 1 



>1o(^,0), 



for R > Zc- The potential produced by the shell satisfies 

V2$,(r, t) = 47rGS(0, t)5[r - R{e, t)]. 
Finally we need the equation of motion for the core mass: 



Zr 



\ dr 



r=Zc,e=0 



Consider linear perturbation modes associated with spherical harmonics Yiq: 

R{e,t) = Ro{t)[l + aiit)Yio], 
V±{e,t) = Ro{t)bi{t)V^Yio, 
S(^,t) = So(t)[l + q(t)llo]. 



(B8) 



(B9) 



(BIO) 



(Bll) 
(B12) 
(B13) 
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Solving equation (p39|), the shell potential to linear order in ai, q is 



^ , , GMs 1 fGM,\ fRo^^+^ 

^ ' ^ r 21 + 1 \ Ro 

for r > R{9,t) (outside the shell), and 

GM, 1 rGM, 



q + (/ + 2)ai 



Y, 



10 J 



$s(r,t) 



r y 



ci- {I- l)ai YiQ 



Rq 2t + 1 \ Rq J \ Rq } 
for r < R{9,t) (inside the shell). Using the unperturbed solution (eqs. 



and 



perturbation equations (p5|)-(B7) for the shell reduce to 



2tai + tci--l{l + l)bi = 0, 



9rdi + 12tai - 2ai 



AMr 



M 



■SllZc + 



AMr 

'm 



2M, 



-ai 



{21 + l)Af 



/(/ - l)ai + -Q 



M, 



Stk + h = ^6nzc + (2^ ^ ^^^^ 



s / -J 



where 6ii is Kronecker delta, and 



4^y/v^c 



i?0 



The equation of motion for the core mass, equation ( BlOj ), becomes 

2M, 



m'^Zc + I2tzc - 2Zc 



3M 



Cl. 



Setting 

ai, bi, Cl, Zc(x {-ty, 

equations ( pi6p -( p20| ) reduce to a set of algebraic equations from which the eigenvalue, 
corresponding eigenmode can be determined. We discuss these eigenmodes below. 



(B14) 

(B15) 

), the 

(B16) 

(B17) 

(B18) 
(B19) 

(B20) 

(B21) 
s, and the 



B.2. I = 1 Modes 

There are six roots for s. Two of these, s = 1/3, —2/3, are trivial modes which do not involve 
any surface density perturbation (ci = 0), and for which the core experiences no acceleration 
{Zc = 0). These correspond to the collapse of a uniform shell onto a displaced core; for 
AR = R- Ro(x i-t), AR/Ro oc (-t)^/^ and for AR = constant, AR/Ro oc (-t)-^/^. 

Two of the remaining four roots correspond to stable modes (s > 0). The two unstable modes 
are shown in Fig. 5. Both lead to the growth of the separation between the center of mass of the 
shell and the position of the core: (i) "Bending" Mode: the surface density perturbation grows 
because one side of the shell collapses faster than the other, and the geometric center of the shell 
moves in opposition to the motion of the central mass; when Mg — > 0, the mode has s = —1 
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and (ai,6i,ci) = (1,0,-2); (ii) "Jeans" Mode: the surface density perturbation grows due to the 
internal tangential flow in the shell, while the shell's geometric center suffers little displacement 
with respect to the position of the central mass; when Mg -^ 0, the mode has s = —1/3 and 
(ai,6i,ci) = (0,1,-4). The "bending" mode grows more rapidly, with s ranging from —1 for 
Ms ^ to s = -(\/T7 + l)/6 = -0.854 for Mc -^ 0. 



B.3. I = 2 Modes 

For / = 2, the core mass experiences no acceleration. There are two types of unstable modes 
as shown in Fig. 6: (i) "Bending" Mode: the north pole and south pole of the shell collapse faster 
and have higher density than the equator, leading to a quadrupolar density perturbation; (ii) 
"Jeans" Mode: the density perturbation is mainly due to tangential fluid motion within the shell. 
The "bending" mode is the more rapidly growing mode, with s in the range between —1 and 
-0.85; 



B.4. Large-/ Limit 

Our results for general / will not be presented here. But it is of interest to consider the I ^ 1 
limit. 

(i) "Bending" Mode: The dispersion relation of bending waves on a pressureless, nonrotating 
surface is o;^ = 27rGEo|/c|, with |A;| ~ l/Ro- This gives 

The wave evolves as exjp{i Jlj dt) = {—t)'. To obtain the amplitude evolution, we use the 
conservation of wave action (energy per unit mass divided by frequency) oc lo\AR\'^, which gives 
\AR\ oc (-t)^/2^ and \AR\/Ro oc {-t)-^/^. Thus 

A /? 1 

-^^{-ty, s = --±isi. (B23) 

/to 

This agrees with our numerical results. 

(ii) "Jeans" Mode: The dispersion relation of density wave is w^ = — 27rGSo|A;|, which gives 
to = zizisi/t. Similar to (i), we find 

A /? f 1 

oc(-t)-l/6exp(i/u;dt)~(-t)^ s = --±si. (B24) 

Ro J 6 

This also agrees with our numerical results. 
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B.5. Effect of Internal Pressure 

We can include the effect of fluid pressure by adding a term — (1/S)V_l(SC|) to the tangential 
Euler equation (p7|); the radial equation ( p6| ) is not affected by pressure. For definiteness, we 
parametrize the shell- averaged sound speed, Cg, by 

Cs = P (^ j , (B25) 

where /3 is a constant. This amounts to adding a term —0^ci to the right-hand-side of equation 

PUD- 



Figures 5 and 6 show the effect of pressure on the eigenvalues of the / = 1 and I = 2 modes, 
respectively. We see that pressure always tends to stablize the modes. However, for Ms/Mc less 
than a few, the "bending" mode is only slightly affected. 

In the large-/ limit, the "bending" mode is unaffected by the pressure, thus equation ( [B23| ) 
still applies. For the "Jeans" mode, the dispersion relation is uj"^ = k'^C'^ — 2TiGTiQ\k\, which gives 

For / ^ Ms/{20^M), we have uj ~ ib(v2/3)/3//t. Using similar precedure as in §B.4, we obtain 

^ « (-t)^ with sc^-\± t^pi. (B27) 

Hq 
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Fig. 1. — Evolution of an / = 2 perturbation during the collapse of a centrally concentrated dust 
sphere. The velocity potential perturbation, 6u, in the upper panel, and the fractional density 
perturbation, Sp/p, in the lower panel, are plotted against the Lagrangian radius, rmit), for three 
different mass shells: m = 0.2 (solid lines), m = 0.6 (short-dashed lines), and m = 0.9 (long-dashed 
lines). The dotted line in the lower panel shows the asymptotic scaling relation as derived in §3.2. 
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Fig. 2. — Evolution of the unperturbed flow during the collapse of a pressure-depleted 7 = 4/3 
polytrope. The adiabatic index is reduced to 7 = 1.3 after the collapse starts. The density and 
velocity are plotted against the Lagrangian radius, Vmit), for of three different mass shells: m = 0.2 
(solid lines), m = 0.6 (short-dashed lines), and m = 0.9 (long-dashed lines). The dotted lines 
show the asymptotic scaling relations. Note that p oc r^'^'^'^~'^' applies to the outer region of the 
self-similar flow (Yahil 1983), while p oc r~^'^ and v oc r^^/'^ apply to the inner region of the 
post-collapse flow. 
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Fig. 3. — Evolution of the / = 2 perturbation during the collapse of a pressure-depleted 7 = 4/3 
polytrope. The unperturbed flow is depicted in Fig. 2. The velocity potential perturbation, 5u, 
in the upper panel, and the fractional density perturbation, Sp/p, in the lower panel, are plotted 
against the Lagrangian radius, rm{t), for three different mass shells: m = 0.2 (solid lines), m = 0.6 
(short-dashed lines) and m = 0.9 (long-dashed lines). The initial perturbation is chosen to be 
5p/p = 1 and 6u = 0. The dotted line shows the asymptotic scaling derived in §3.2. 
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Fig. 4. — Same as Fig. 3, except that the initial perturbation is chosen to correspond to the gi 
mode of a 7 = 4/3 polytrope with adiabatic index 71 = 5/3. 
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Fig. 5. — Perturbation modes (/ = 1) in a collapsing shell. The solid curves are for zero pressure 
(/3 = 0); the lower curve corresponds to the "bending" mode, and the upper to the "Jeans" mode. 
The dotted and dashed curves include pressure -with (3 = 0.5 and /3 = \/0.5, respectively. Note that 
when s is complex, only its real part is plotted. 
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Fig. 6. — Perturbation modes {I = 2) in a collapsing shell. The solid curves are for zero pressure 
(Ji = 0); the lower curve corresponds to the "bending" mode, and the upper corresponds to the 
"Jeans" mode. The dashed curves include pressure with /3 = VO.S. Note that except for the 
lower-left branch, the dashed curves depict only the real part of s. 



